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The physics of hadrons, studied with electron scattering, is an active research field: Jefferson Lab is completing 
its successful program using a 6-GeV CEBAF accelerator and is preparing a 12-GeV upgrade to start a new physics 
program within the next several years Electron-scattering experiments probing the structure of hadrons also 
continue at the MAMI and ELSA facilities in Germany. A common feature of the modern experiments that study 
electromagnetic structure of hadrons is their high precision. 

The increased accuracy of the measurements put stringent requirements on QED radiative corrections that must be 
applied to the data. The largest QED corrections arise from bremsstrahlung, i.e. the emission of real photons by the 
electrons. The bremsstrahlung corrections are characterized by large logarithms of the type logQ 2 /ra 2 , where Q is 
transferred 4- momentum and m e is electron mass, resulting in an order-of-magnitude enhancement of 0{a em ) effects. 
Systematic uncertainties due to the bremsstrahlung corrections can be minimized by using iterative procedures and 
without requiring additional knowledge of hadronic structure. However, when the experiments reach accuracy at 
a percent level, their interpretation also becomes sensitive to box-type contributions described by two-photon (27) 
exchange. A good example of the importance of 27-exchange corrections is the Rosenbluth separation experiments 
in elastic electron-proton scattering [2|. In these experiments, the measured Rosenbluth slope may be misinterpreted 
due to the omission of the angular dependence of 27 corrections in the data analysis, which can mimic a significant 
part of the electric form factor contribution to the cross section. 

In this paper, we evaluate 27-exchange effects for the processes of exclusive pion electroproduction e+p — > e+n+n + 
and e + p — > e + p + ir°. These processes have been studied experimentally in great detail and have become a valuable 
source of information on the electromagnetic structure of nucleons. For example, depending on the choice of kinematics, 
they provide access to the electromagnetic transition amplitudes from the nucleon ground state to excited states @; 
or Q 2 -dependence of the pion electromagnetic form factor |^-|7|; or, in a deep- virtual regime, polarization-dependent 
Generalized Parton Distributions of a nucleon Q . 

The QED radiative corrections are indispensable in the interpretation of the experimental data on electron scatter- 
ing. Although significant theoretical effort has been dedicated to this problem, more work still needs to be done. A 
classical approach developed by four decades ago was mainly intended for inclusive and elastic electron scattering, 
assuming integration over the phase space of final hadrons. More recent work addressed the case of electroproduction 
with detection of hadrons in the final state for semi-inclusive [ll[ and exclusive diffractive [l2j processes. A calcu- 
lation of QED corrections to exclusive electroproduction of pions by [13] provided the next-to-leading order QED 
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radiative corrections to the cross section and the beam polarization asymmetry for exclusive pion electroproduction. 
In comparison to earlier work , the authors of Ref.[13j (a) eliminated dependence on a nonphysical parameter 
coming from splitting soft and hard regions of the phase space of radiated photons by covariant treatment of an 
infrared divergence using the method of [l4[, O 3 ) treated bremsstrahlung exactly, without applying either the soft 
photon or peaking approximations and (c) retained the effects arising from the kinematics of detection in coincidence 
of two particles, an electron and a hadron, in the final state. However, the calculation of |13l | only included vacuum 
polarization corrections and contributions from bremsstrahlung and vertex corrections from the electron. No box 
diagrams with two-photon exchange or photon emission from hadrons were considered. 

Motivated by the need to evaluate two-photon exchange for measurements of the pion form factor, the authors of Ref. 
[l5j extended their approach to elastic electron- proton scattering (outlined in [l6j|) to calculate two-photon exchange 
corrections in the elastic electron-pion scattering, including only pion elastic intermediate states. The authors of 
Ref. [l7| considered a process of eir + elastic scattering as well. Both these papers, pj| and [It]], concluded that the 
two-photon exchange plays a sizable role in the pion form factor measurements at extreme backward angles. In [l8j . 
the two-photon exchange amplitude for the elastic electron-pion scattering was computed in the dispersion approach 
and with the monopole parameterization of the pion form factor, including both elastic and inelastic contributions. 
All the previous calculations of 27-exchange for pion form factor measurements treated the pion as a stable and free 
particle at rest, which complicates the application of this correction to actual measurements. 

In the work presented here, we calculate the two-photon box corrections for the cross sections of the entire process 
of pion electroproduction on a nucleon for two production channels, 

e(fci) + p{k 2 ) e(fc 3 ) +7T+ (fc 4 ) + n{k 5 ) 
e(k 1 )+p(k 2 ) -> e(k 3 )+p(k 4 )+7T°(k 5 ) 

where in the final state an electron and a charged hadron are detected, 7r + or the proton. The calculations are per- 
formed both analytically and numerically using FeynArts and FormCalc [19| as the base languages. The regularization 
of infrared divergences arising in the box diagrams is addressed by giving the photon a small rest mass and canceling 
that by adding the soft-photon bremsstrahlung contribution from the electron, proton, and the charged pion. 

The paper is constructed as follows: In Section II, we derive formulae for 27 corrections to pion electroproduction 
amplitudes. 

Section III describes the soft-photon bremsstrahlung contribution required for cancellation of infrared divergences. 
Section IV provides the analytical results. Numerical analysis is performed in Section V, and the conclusions are given 
in Section V. 



II. TWO PHOTON BOX DIAGRAMS IN PION ELECTROPRODUCTION 

Our goal is to derive the radiative correction associated with the two-photon box diagrams in exclusive pion 
electroproduction. Here, we choose to employ, compare and analyze two approaches. The first approach is based on 
the soft-photon approximation, where we follow two prescriptions. In the first prescription, we apply the soft-photon 
approximation in the box amplitude consistently in the numerator and denominator algebra according to the method 
suggested by Tsai [20]. In this paper, we will call this type of soft-photon treatment of the box amplitude SPT 
(Soft-Photon- Tsai). The second prescription is based on the work of Maximon and Tjon [21|, where the soft-photon 
approximation is applied only in the numerator algebra. Accordingly, we will refer to this prescription as SPMT 
(Soft-Photon-Maximon-Tjon). In the second approach, we evaluate the box amplitude exactly, without relying on 
the soft-photon approximation. Since in this case the box amplitude can not be factorized by the Born amplitude, 
this approach will result in the model dependence of the radiative correction. Here, we evaluate the box amplitude 
by employing the model with the monopole form-factor in the photon-proton, photon-pion and pion photoproduction 
couplings. We will refer to this approach as FM (Formfactor- Model) . 

The two-photon box diagrams in exclusive pion electroproduction are given in Fig.([T]). For the incoming electron 
and proton, we have momenta k\ and k 2 . For the outgoing electron, detected proton (or pion: tt + ) and undetected 
neutron (or pion: ir°) we have momenta k 3 , k 4 and k$, respectively. 

Let us start with the definitions of the input kinematic parameters. We choose to use the incoming electron energy 
in the laboratory frame Ei a t,, the momentum transfer Q 2 = —q 2 = — (£3 — ki) 2 and the invariant mass of the virtual 
photon plus incoming proton W 2 — (q + k 2 ) as input parameters. In addition, the polar (#4) and azimuthal (^4) 
angles of the detected proton (pion) are defined in the center of mass of the final hadrons. The extended set of Lorentz 
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Figure 1: The two-photon box diagrams in pion electroproduction. 

invariants is defined according to the prescription of fl3^ . 

S = 2 (kik 2 ) = 2m 2 Ei ab , 

X = 2 (k 3 k 2 ) = S - W 2 - Q 2 + ml, 

S t = 2{k 2 k i )=W 2 + m 2 -m 2 + V 3 -V ll (1) 
Vi t 3 = 2(ki^ 3 k i ) = 2(E 1 ^E i - pi^Pi (cos 64, cos #1,3 + sin 64 sin #1,3 cos ^4)), 
Ul = S-Q 2 ,u 3 = X + Q 2 , 

Ai,3 - u? i3 - 4m?T^ 2 , A 9 = (S - X) 2 + Am 2 Q 2 . 

With the help of the invariants in Eq.([T]), we define energy and momenta in the center of mass system of the virtual 
photon q = k\ — k 3 and initial proton as follows: 

77 _ u i,3 _ \Al,3 



S - X + 2m\ J\ 

E 2 = — > P2=Pq = 



2W 



2W 



(2) 



E 4 = 



W 2 +m(- 
2~W~ 



-, Pi 



J (W 2 + m\ - ml) 2 - AmlW 2 



2W 
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The angles #1,3 are given by the following expressions: 



a 



1,3 (S - X - 2Q 2 ) ± 2Q 2 W 
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00801,3 = — ! 7= — 7= , (3) 



2WV\ , , ., 

sin 6*1,3 = 7= , X = Qu 1 u 2 -QW-mfX q 



A. Soft-photon Approximation 

The soft-photon approximation significantly simplifies the calculations of the two-photon box amplitudes. Specifi- 
cally, a photon which does not couple to the pion photoproduction vertex is treated as a soft one. This implies that 
we extract only the infrared-divergent part of the amplitude and, for now, we assume that the non-infrared-divergent 
part of the amplitude gives a negligible contribution. The infrared divergences arising in the cross section from the 
two-photon box diagrams will be treated later using soft-photon bremsstrahlung. 

As an example, let us start with the first graph on Fig.([3Ja)). For this graph, we can write the following: 



Mi = 



u(k 3 ,m 3 ) (iej^) — — l - (iej„) u(fci,mi) 



16tt 4 

(fc 5 (4),m 5(4) )(ir 7 -^" (p) ) 



(fci - q) - m\ 

/&2+ A + m 2 



(k 2 + q) -mi 



(r^- p (q))n(k 2 ,m 2 ) 
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g2 ' ( g + fc 3 _ fcl )2- 



(4) 



Here, = m n for the case of 7r + production and 7714 = m p for tt° production. The photon-proton coupling 
has the following structure: r P p~ 7 ~~ p (q) = —ie^fpF(q), with the form-factor F(q) = A ^_ q 'i (A 2 = 0.83m 2 ,) taken 
in the monopole form. In the case of the soft-photon approximation with either the SPT or SPMT (q — > GeV) 
prescriptions, the coupling T P j ~ 1 ~ p (q) converges to the coupling of the point-like proton r^~ 7 p (q — > 0) = —iejp, and 
the box amplitude becomes independent of the choice of the type of form-factor. We do not show the explicit structure 
of the pion photoproduction coupling Ta 7 n here because the structure of that coupling has no impact on the 
two-photon box radiative correction calculated in the soft-photon approximation. 

Using the soft-photon approximation according to the SPT prescription, we neglect the momentum of the photon 
in the numerator algebra and in the denominator of photon propagator: 

help of the Dirac equation, terms such as (jk\ +m t ) (ier/ v ) u(ki, rrii) can be simplified into 2ie(ki) u - For the amplitude 
in Eq.((4]) we can write: 



, , ? — : — rr — > ., 9 , . As a result, with the 

(8+K3— *i) (K3-K1) 



Ml 



SPT 



[u(k 3 ,m 3 ) (ie7 M ) u(fci,mi)] w(fc 5 ( 4 ), m 5 ( 4) ) \T p a 7 w u(k 2 , m 2 ) 

4(fcl-fe 2 ) j 



fj 



ft-CX 



(h-hy 



(k 2 +q) -m 2 (ki - q) 
Taking into account that the Born-level amplitude is 



M Q = [u(k 3 ,m 3 ) (ze7 M )u(fci,77ii)] u(fc 5 ( 4 ) , m 5 (4)) (j** 7 * " (p) ) u(k 2 , m 2 ) 



(fc 3 -fci) 2 



we can write 



(5) 



(G) 



SPT 



2tt 



S ■ Cq ({fci,mi}, {-k 2 ,m 2 }) ■ M . 



Here, Co ({ki, m,}, {kj, rrij}) is the Passarino-Veltman three-point scalar integral defined as: 



C ({ki,mi},{kj,mj}) = -r-? f d 4 q-^ 

i7r J Q 



(ki - q) - mf (kj - q) - to 2 



(7) 



(8) 
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M 2 SPT = 


a 






M 3 SPT = 


a 






Ml PT = 


a 




~2n 



A similar approach can be employed for the rest of the graphs in Fig.JT]). The results for graphs of types (b), (c), and 
(d) can be written in the following form: 

" X ■ C {{k 3 , m 3 }, {k 2 , m 2 }) • M , (9) 
V 3 ■ C {{k 3 ,m 3 },{-k 4 ,m 4 }) ■ M , (10) 
Ki-Co({fci,mi},{fc4,m 4 })-Afo. (11) 

Combining Eqs.Q, ®, (fT0|) and (fTT|) we can write the total amplitude for the two-photon box diagrams in the 
following form 

4 

Mt P J = Y, M i =-^ { - S - C ° m i}' {" fc 2, rn 2 }) +X-C ({k 3 , m 3 ], {fc 2 , m 2 }) + 

^•Co({fc3,m 3 } ) {-fc4,m4}) + ^-Co({fci ) mi},{fc4,m4}))-Mo. (12) 

It is worth mentioning that in the case where ir + appears in the loop (FigJT] (c, d)) of the box diagram, we use the 
7r — 7 — 7r coupling in the following form: = — ie ■ (fc™ + k° ut ) a F(q). The pion propagator is given by 11^ = 

b2 l m i ■ When the soft-photon approximation is applied, the photon-pion coupling becomes _^ —2ie ■ (k 4 ) a 

and hence the box amplitude will have the same general structure as for the case of ir° electroproduction. 

It is obvious that in Eqs.([7]) and (|9"l lll[) , the box amplitude is factorized by the Born amplitude, which can be 
evaluated in the different models or approaches. Keeping terms of order 0(a 3 ) only, the factorization of the pion- 
electroproduction cross section can be accomplished in the following way: 



da 



SPT 



(|M | 2 + 2Re [Ml P jMl\ + O (a 4 )) dT 2 ^ 3 = da Q + da s h p x T = da • (1 + 8§^) , (13) 



SPT 

box ~ dao 



where dao is the Born differential cross section, dT 2 ^. 3 is the phase factor for the 2 — > 3 process and Sf PT — (7hox 
is the radiative correction due to the box diagrams contribution: 

8bo P J = -~Re[S ■ C {{h, mi}, {-fa, m 2 }) + X ■ C ({k 3 , m 3 }, {fc 2 , m 2 }) + 



2] 

7T 

V 3 ■ Co ({fc 3 ,m 3 },{-fc 4 ,m 4 }) + Vi ■ Co ({fci,mi},{fc 4 ,TO4})]. (14) 



Since the dependence on the pion photoproduction coupling rP~T _7r_Tl (p) j s cance i ec j j n the soft two-photon box 
radiative correction, we can call this type of calculation a model-independent type. 

The three-point scalar integral Co ({fcj, m^}, {kj, rrij}) in Eq. (|14[) was evaluated both using an approximated ex- 
pression (nii -C 771 j ) taken from [22], and exactly with the help of the package LoopTools [l^]. The explicit expression 
for this integral reads as follows: 



C ({ki,mi},{kj,mj}) = - 



4(fei • kj) 



In \ J ' + 2 In v 2 J; +ln^ In 

rrij \ uij m i I * 



l ln2 m|_ _ m j-2(k l -k J ) 
2 m? I 2(h-kj) 



(15) 



where Li 2 (z) = — j Q z ln ( 1 ^~ t ) fa j s the usual dilogarithm function. It is important to note that in our calculations we 
did not apply the approximation Re [Co ({A;*, mj}, {— kj, rrij})] ~ —Co {{ki,mi},{kj,rrij}) introduced earlier in [20]. 
As a result, this produced a shift of an in the value of the box radiative correction using the SPT prescription. 

When evaluating the two-photon box amplitude using the SPMT prescription of the soft-photon approximation, 
we follow the same steps as in SPT, but we no longer neglect the photon momentum (q) in the denominator of the 
photon propagator ^ + ^_ fc ^ 2 . This effectively results in the replacement of the three-point scalar integrals by the 
four-point ones. Specifically, we can write: 

SZ MT = -—Re[S ■ Do (fo, mi}, {-fa, m 2 }, {fa,k 3 }) + X ■ D Q ({k 3 , m 3 }, {k 2 , m 2 }, {-fa, -k 3 }) + 

TT 

+V 3 ■ D ({k 3 ,m 3 }, {-k 4 ,m 4 }, {-ki,-k 3 }) + Vi ■ D ({fa, mi}, {fa, m 4 }, {fa, k 3 })] . (16) 
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The four-point scalar integral Dq has the following structure: 

D ( {kl , m lh {kj ,m j} , {k a ,k,}) = ■ (fci _ 9 j 2 _ m? • (fe ._^_ m? ' + 

and can be numerically evaluated using the package LoopTools. To see how our results in Eq. (|17l) compare to the 
results of [2l|, we have reduced Eq. (fT()]) to the case of elastic electron-proton scattering and for that the corresponding 
correction can be written in the form: 

SZ-(e P) = ~—Re[S ■ D ({fci.mi}, {-fc 2 , m 2 }, {h, k 3 }) + X ■ D Q ({fc 3 , m 3 }, {k 2l m 2 }, {-fci, -fc 3 }) ] • (18) 



Here, since this is elastic e—p scattering, we take — (GeV) and invariant mass W = m p , which effectively reduces 
the kinematics of the 2 — > 3 process to those of a 2 — > 2 process. After the treatment of infrared divergences in Eq. (fT5)) 
with soft-photon bremsstrahlung, we find that the results produced in Eq.(ll) from (2ll | (fo r the case of 7 — 7 box only) 
and Eq.(|18| are identical. Such an approach was considered by Maximon and Tjon [21] as a "less drastic" version 
of the soft-photon approximation. We note, however, that since it treats numerators and denominators of fermion 
propagators differently it violates Ward identity, i.e. it does not preserve the gauge invariance. In principle, such an 
approximation would lead to the non-renormalizability of the theory. Luckily for the specific case of the two-photon 
exchange calculations, it does not result in any unphysical divergences. 



B. Exact Model-dependent Approach 



In the second approach, we apply the computational methods described in (l9| and calculate the two-photon 
box correction without relying on the soft-photon approximation. In this case, the two-photon box amplitude is 
calculated exactly, but since it is not possible to factorize it by Born amplitude, this approach becomes effectively 
model-dependent. Here, we define a radiative correction using the Formfactor-Model (FM) notation in the following 
form: 



j„fm 2Re 

cFM _ aa box _ 
°box — 



da I Mo I 



(19) 



In order to calculate the radiative correction 8™J to the Born cross section, we employ a model where the monopole 
form-factor is applied in every photon-hadron coupling. In this case, the pion photoproduction coupling is described 
by the axial- vector contact term 7 71 n — e ^ +F ^ 7 a 75-F(g), which is weighted by the form-factor taken in 
monopole form F(q) — ■ The low-energy constants D = 0.40 ±0.03 and F = 0.61 ±0.04 were determined earlier 

in [23| and — \[2 ■ 92 MeV is the pion coupling constant. The photon-proton and photon-pion couplings also have 
the same structure as in the case of SPT or SPMT: r^ 7_p (q) = -iejpFiq) and rj-T-^ = _j e . + k° ut ) a F(q), 
respectively. We do not neglect the momentum of the photon in the form-factor when evaluating the loop integral. 
Calculations of the two-photon box radiative correction where done analytically using FeynArts and FormCalc and 
later numerically using LoopTools. Although we have analytic expressions derived for the correction in Eq. (|19|) as 
well, it is not possible to show them in this paper due to their length. For this approach, we show numerical analysis 
only, which is presented in Section V. 



III. SOFT-PHOTON BREMSSTRAHLUNG 



At this point we are ready to address the infrared divergences arising in the box diagrams. Regularization of the 
infrared divergent integrals in Eq.© and (fl~7j) can be done by assigning a small rest mass A to the photon. Clearly 
the photon mass should be removed to avoid nonphysical representation of the final results. This can be achieved by 
adding the soft-photon bremsstrahlung contribution coming from the graphs in Fig.Q. 

As an example, we consider the first graph in Fig.© and evaluate the amplitude in the soft-photon emission 
approximation by neglecting the momentum of the photon in the numerator algebra. Hence, we can write: 



u(k 5 ,m 5 ) ( 



pP-7" 



Mi i7 = i [u (k3,m 3 ) (iej^) u(fci,mi)] 



(20) 
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n(7c°) 





Figure 2: Representative graphs for the bremsstrahlung process in pion electroproduction. Here, the momenta distribution is 
the same as in the case of box calculations with momenta of the photon taken as fco. 



Here, e*' a (fco) is the polarization vector of the emitted photon. In the limit where the energy of the emitted photon 
is small, we can readily assume that rP~ 7_p (fco) — > rP _7_p (0) = iQ p ^ a and (^Zk'^^m 2 (* ( 9p7a) w (^2, "12) can be 

(fca)c 



02-fco) 

reduced to —iQ p ^^"-. Mfe; ™ 2 ) with the help of the Dirac equation. As a result, we can write the following for the 



amplitude in Eg. (1151) : 



Mi, 7 = M ■ Q P 



(k 2 e* (fco)) 



(fc 2 ■ fco) 

We can write bremsstrahlung amplitudes for the rest of the graphs in a similar way: 



Mi, 7 = M 



. foe* fa)) 
(fci • fco) 



(fc 2 e* (fc )) _ (fc 3 e* (fc )) 



(fc2 ' fco) 



(fc3 ' fco) 



(fc 4 £* (fcp)) 

' T+(p) (fc 4 -fco) 



(21) 



It is obvious from the equation above that the cross section for the bremsstrahlung process is parameterized by the 
Born cross section as: 



da*. 



da 



1 



(2tt) 3 



d 3 k 
2tu 



, (fcie*(fco)) 
(fci • fco) 



Qp 



(fc 2 e* (fc )) (fc 3 e* (fco)) 



(fc2 ' fco) 



{h ■ fco) 



Qt 



(fc 4 e* (fco)) 



(fc 4 • fc ) 



(22) 



Here, lu is the energy of the emitted photon. For the calculations of the bremsstrahlung cross section, we have to 
account only for the parts of Eg. (1221) which are responsible for the treatment of infrared divergences in boxes. Namely, 
the products of graphs one and three, one and four, and two and three plus two and four in Fig. ([2]). After summing 
over all polarization states of the emitted photon and using (fc^e (fco)) (fcj e * (fco)) = — (fcjfcj), we can write, for the 
part of the bremsstrahlung cross section responsible for the IR-divergence cancellation in boxes, the following result: 



a 

'2^2 



d 3 kn 



(fci • fc 2 ) 



dcr 7 = da ■ Ay, 

(fci • fc-l) 



(23) 



(fc 3 • fc 2 ) 



(fc 3 • fc 4 ) 



(fci ■ fc ) (fc 2 ■ fco) (fci • fco) (ki ■ fco) (h ■ fc ) (fc 2 • fco) (k 3 ■ fc ) (fc 4 ■ k ) 



where <5 7 is dehned as the soft-photon bremsstrahlung correction. Here we have used Q w +(p) = — e. The integral in 
Eq. (|23p clearly has IR-divergent behavior ,and its regularization can be achieved in the same way as before by giving 
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the photon a fictitious mass (A) with the photon's energy being defined as u> = ^/k 2 , + A 2 . Eq. (f2"3"| is presented in an 
invariant form, so integration over the phase space of the emitted photon can be carried out in any reference frame. 
We choose to follow the prescription of [2l| where integration is carried out in R frame defined as ks = — k . With 
the final results written in invariant form, we list only the final equations relevant to this work. The soft-photon 
integral can be written as follows: 



d 3 k Q 



1 



w (fci • k ) (kj ■ ko) 



The integral I (ki,kj) is evaluated as: 



where for i ^ j we get 



I {ki-i kj 



2tt 

Hj 



A 



(i) 



A 



(2) 



A w = i n ( a? ,^V h / 4Ae 



m, 



and 



.(2) 



In' 



miVA 2 / 



-Li-, 1 



uij v A 2 



Li 2 [l 



V A 2 

A (kj ■ A) 
A 2 7» 



-Li 2 \l- 



a ijPjlij 



where 



Otij ki hj 



A = (ki ■ A) + yj(ki- A) 2 -m 2 A 2 ; 7ij = ^ (h ■ k 3 ) 
In the case where i = j, the integral in Eq. (|24p has the simple form: 



(24) 



(25) 



(26) 



Li 2 1 



(jjj ■ A ) 



(27) 



2 2 2 

infm 1 * 
L j 



1 (kifki) y 



hi 



/4Ae 



V A 2 



2 (fcj • A) 



In 



2 A 2 



(28) 
(29) 

(30) 



The momentum A is determined by using the soft-photon approximation and is equal to the momentum of undetected 
particle A = k§. The Ae parameter is the upper value of the energy of the emitted soft photon. If hard-photon 
bremsstrahlung is considered, the dependence on Ae is canceled and replaced by the experimental cuts. In this work, 
we have decided to keep only the soft-photon contribution and determine Ae from the pion production threshold 
conditions. That translates into the requirement that the invariant mass of the undetected hadron and the emitted 
soft photon should not exceed the value of + m^. As a result, the maximum energy of the emitted soft photon in 
R frame can be calculated in the following way: 



Ae 




where x — —. 

m 5 



(31) 



For the case of ir° electroproduction, we have = m n , and therefore Ae = \m n . For 7r + electroproduction, we have 
m 5 = m n and with the help of the approximation x <C 1 we can write Ae = m v . 

The scalar products relevant to the soft-photon integral are written in terms of the invariants of [24] and are defined 
as: 



(fci • k 5 ) = Ul 1 , (k 2 ■ k 5 ) 



(&3 ' h) 



2 

U 3 - V 3 



{ki ■ fc 5 ) 



S-St-X 



V!-V 3 + s t 



(32) 



9 



Finally, we can express the soft-photon factor in terms of the evaluated integral I (fej, kj): 



S ~t = -^2-(- S - 1 ( fc i> fc2 ) + Vl ■ 1 ki)+X-I (k 3 , k 2 ) -V 3 -I (h, fc 4 )). 



(33) 



A combination of the box corrections for any of the cases in Eqs. (fH|) . (11611 and (fTO|) and the soft-photon bremsstrahlung 
factor in Eq. (|33p cancels the photon mass analytically and makes the total correction free of this nonphysical parameter. 



IV. RESULTS 

In this section we provide the analytical results for the box correction by deriving compact equations using the SPT 
prescription only, and show the explicit cancellation of the photon mass A. For the cases of SPMT and FM, analytical 
expressions for the four-point tensor integrals are quite lengthy, so for these cases we will only show numerical results. 

Let us define the total SPT differential cross section and the box correction treated with the bremsstrahlung 
contribution as: 



dag™ = da SPT + m = d(Jo (1 + xspt + ^ = d(Jo (1 + 6 sp T) 



(34) 



and separate the infrared-finite part Sf ot from the part regularized by the photon mass 6*$ in the following way: 

6*F = C + 8[ ot . (35) 



Here, we can write = + <5^ R , where the infrared-regularized corrections 8^. and 8* R have the following simple 



structure: 



X™ - --In ^ 



X V 3 



(36) 



a, 4Ae 2 

7T A 1 



1 S ^1 Vl 

ln x +ln v 3 



(37) 



It is clear now that a combination of Eq.(|31)|) and Eq. (|37p results in the finite total correction which can be written 
in the following way: 



xiR _ xiR , xi R - i n m 2 

°tot ~ °box + °7 - n 111 4Ae 2 



In In — 

X V 3 



(38) 



The infrared-finite part of the correction 5f ot = 5^ ox + 8^ can be simplified for the box part as 



8 F 

°box 



a 

7T 



1 , S , S ■ X 1 , V 3 , Vi ■ v 3 



In — - In 

2 X 



In — • In ■ 

2 Vi 



7T — Lio 



S + mi 



Li 2 



X-m\ 
X 

V 3 



(39) 



and the soft-photon finite bremsstrahlung part has the following structure: 



/Vfr-viA + Ll2 L _ ml-^-v^ _ U2 (l _fk. ( Ul - Vi) 



S-ml 



Li 2 1 



m| ■ (m - Vj) 
Vi-/3 4 



T . , . #2 • (U3 ~ Vj) '. , - . , , 

Ll2[1 - X-m\ 1 £ ' jl ' 



Lia 1 - 



A ■ («3 ~ V 3 ) 

V3 • ml 



Lia 1 - 



Vi • m§ 

TO^ ■ (tt 3 - V 3 ) 
X-p 2 

mj ■ (u 3 - V 3 ) 
V 3 -/3i 



(40) 



10 



Here, functions 




j8 4 = m 4 + 4/ mf — mf • mf . (41) 



A combination of equations Eqs. (|35l) . (f3T?)) and (|4"0")) represent a rather compact set of final expressions for the 

box correction to the pion electroproduction evaluated using the SPT prescription of the soft-photon approximation. 
At this point we have everything ready to proceed to the numerical analysis. 



V. NUMERICAL ANALYSIS 
A. The 7T° electroproduction 

Let us start with the numerical implementation of the analytical expressions derived for the two-photon box radiative 
correction for the kinematics relevant to the recently completed CLAS [25] experiment performed in Hall C at JLab. 
The goal of this experiment (25| is to measure the N — A transition form factors in the region of high momentum 
transfers. Previously, it was found that the radiative correction arising from the electron vertex corrections, boson 
self-energies plus emission of soft and hard photons from an electron could reach as much as —23% [l^ . Evidently, 
these corrections play a vital role in the extraction of the Born cross section from the measurements. In addition, 
the size of the two-photon box radiative correction could also become an important part of that analysis. Fig. ([3]) 
shows our results for the box correction in the region of A resonance as a dependence on the azimuthal and polar 
angles in the hadronic center of mass (hadronic kinematic variables) for all approaches considered above and for the 
experimental runs with Q 2 = 0.4, 6.36 GeV 2 and electron beam energy Ei a b = 5.75 GeV in the lab reference frame. 
We use the blue dot-dashed line to show the correction calculated using the SPT prescription. The gray dotted curve 
represents the correction obtained in the SPT approach but with subtracted air term. That is equivalent to the 
result produced by the approximation Re [Co ({fcj, mi}, {~kj, nij})] ~ —Co {{ki, mi}, {kj, m ■,"}■), which was introduced 
earlier in [20j . The green dashed curve is a correction obtained using the SPMT method (see (2lj). The red solid line 
corresponds to the result produced from the exact but mo del- dependent FM approach. 

For high momentum transfers the correction in Fig. (top row) shows the strong angular dependencies in the 
SPMT and FM approaches and ranges from —10% to —20%. As for the SPT approach, the correction is almost 
constant and lies in the interval of —0.5% ~ —3.0%, depending on whether we use the approximation of [20] in our 
calculations. In general, from Fig. ([3]) we can see that for high momentum transfer the correction calculated by the 
SPT method is notably smaller compared to the SPMT and FM cases. The absolute difference is large and is in the 
interval of 13% ~ 18% for both the polar and azimuthal dependencies. This discrepancy can be explained by the fact 
that in the SPMT and FM cases the correction is calculated using the base of four-point tensor integrals and therefore 
it is directly enhanced by the scale of the momentum transfer oc — Q 2 In (see Eqll6p. Hence, the difference between 
the SPT and SPMT, FM approaches is substantially enhanced at higher momentum transfers. More specifically, for 
Q 2 = 0.4 GeV 2 , Eiat, — 5.75 GeV and W = 1.232 GeV, we can see that the difference between all approaches is smaller 
and is of the order of 3%. This primarily happens due to the fact that the absolute value of the correction calculated 
using the SPMT and FM approaches is substantially decreased for small momentum transfer (see Fig. ([3]), bottom 
row). With the increase of the invariant mass, we also observe that the correction tends to become larger at higher 
momentum transfers. 

Fig. ((3]), relevant to the earlier CLAS experiments (26l - [2^ |. shows the low momentum transfer Q 2 = OA GeV 2 
two-photon box correction dependencies on the invariant mass and the hadronic polar and azimuthal angles. The 
dependence of the correction on the invariant mass at Q 2 — 0.4 GeV 2 is rather weak for all cases with the exception 
of the SPMT approach. The angular dependencies show that the correction is rather small and lies in the interval of 
-3.0% ~ 2.5%. 

The study of the radiative correction with respect to the leptonic kinematic variables can be reproduced by ex- 
pressing the correction as a function of the virtual photon degree of polarization parameter e, which can be evaluated 
using the following expression: 
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Figure 3: tv° electro-production two-photon box correction angular dependencies for the high Q 2 = 6.36 GeV 2 (top row) and 
low Q 2 = 0.4 GeV 2 (bottom row) momentum transfers, W = 1.232 GeV and Ei a b = 5.75 GeV. Left column: dependence on 
cos #4 with (j>4 = 180°. Right column: dependence on <^4 with 64 = 90°. Dot-dashed curve - SPT, dotted curve - SPT with an 
subtracted, dashed curve - SPMT, solid curve - FM approach. 



Here, 9i a b is the electron scattering angle in the laboratory reference frame, which will replace Ei a b as an input 
kinematic parameter. In this case, E\ a b is no longer an independent variable and can be calculated in the following 
way: 



E lab = — !— \Q 2 + W 2 -m 2 2 + J (Q 2 + W 2 - ml) 2 + ) . (43 ) 

4m 2 ^ y snr ei ab J 

For the fixed high momentum transfers of 3.0 GeV 2 and 7.0GeV 2 , the dependence of the two-photon box correction 
on e is shown in Fig. ([5]) (left and middle plots). Evidently, the box correction calculated in the SPT approach shows 
a relatively weak dependence on the degree of polarization parameter for high momentum transfers and tends to grow 
for the case of backward scattering. As expected, in the SPMT and FM prescriptions, the correction has a strong 
dependence on e and can reach values of —17% ~ —50% for backward electron scattering at Q 2 = 3.0 ~ 7.0 GeV 2 . It is 
interesting to observe that for the high momentum transfer case, the correction calculated in the region of A resonance 
using the FM approach has the strongest dependence on e and the largest value compared to the results obtained with 
the SPT and SPMT prescriptions. This effect could be explained by the fact that for the high momentum transfers 
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Figure 4: it electroproduction box correction for the fixed Q 2 — OAGeV 2 and Ei a b = 1.645 GeV. Left plot: dependence on W 
with <j>4 = 90° and 64 — 90°. Middle plot: dependence on cos #4 with <f>4 = 90° and W = 1.232 GeV . Right plot: dependence 
on (f>4 with 64 = 90° and W = 1.232 GeV . Dot-dashed curve - SPT, dotted curve - SPT with an subtracted, dashed curve - 
SPMT, solid curve - FM approach. 
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Figure 5: 7r° electroproduction two-photon box correction (for detected proton) dependencies on virtual photon degree of 
polarization parameter e for momentum transfers Q 2 — 3.0 GeV 2 (left plot), Q 2 — 7.0 GeV 2 (middle plot) and Q 2 = OAGeV 2 
(right plot). All plots are given for (j> 4 = 90° and 64 = 90° and W = 1.232 GeV. Dot-dashed curve - SPT, dotted curve - SPT 
with air subtracted, dashed curve - SPMT, solid curve - FM approach. 



and invariant masses in the region of A resonance, the two-photon box amplitude evaluated in the FM approach 
becomes insensitive to the monopole form factor. Since we have to divide the interference term in the numerator of 

the box correction (S^f = — ^j^p— by |-^o| 2 > this effectively produces an overall enhancement of the correction 

by a factor of A ~f? ■ I n order to see if that is in fact true, we evaluate the two-photon box correction exactly but 
without the monopole form factor in the couplings, hence treating the nucleon and the pion as point-like particles (see 
Fig.®). A reduction in the correction value due to the point-like nucleon and pion would point out that the two- 
photon box amplitude calculated in the FM approach is not sufficiently suppressed by the presence of the monopole 
fromfactor in order to overcome the — j^i enhancement coming from the denominator of the box correction. 

As can be seen from Fig.©, the correction calculated using the exact approach with the point-like nucleon and 
pion is reduced substantially (it is now even smaller than the correction obtained with the SPMT approach), hence 
justifying the behavior of the correction calculated in the FM approach at the high momentum transfers. For the low 
momentum transfer Q 2 = OAGeV 2 (see Fig.© (right plot)), the correction becomes smaller and all three approaches 
produce similar results in the range of —5.0% ~ 2.0%. 
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Figure 6: 7r°electroproduction two-photon box correction dependencies on the virtual photon degree of polarization e for the 
case of high momentum transfers Q 2 = 3.0 GeV 2 (left plot), Q 2 — 7.0 GeV 2 (right plot) and without the monopole form factor 
in the case of exact calculations (solid line). All plots are given for cj>4 = 90° and #4 = 90° and W = 1.232 GeV . 



B. 7r electroproduction 

The exclusive ir + electroproduction process proved to be a very successful tool in studies of the transition form 
factors to the states in the mass region above A(1232) resonance (see [3j). Since there is a cluster of three states 
#(1440), #(1520), #(1535) and at least nine N* , A* states available in the mass region 1.44 GeV ~ 1.72 GeV, it is 
crucial to distinguish them in order to be able to gain information about the transition form factors. Most of these 
states are quite sensitive to both mr + and pn channels, but since the states with isospin 1/2 couple more strongly 
to the n7r + than to the pn , it is very important to obtain a full picture about the 7r + electroproduction cross section 
and extract the unradiated cross section. Obviously, a full picture requires the inclusion of the radiative corrections. 
Here, we have extended calculations of [13j with the two-photon box correction for the ir + electroproduction case. As 
an example, in the Fig. (JTJ) (left and middle plots), we show the box correction dependencies on the photon degree of 
polarization parameter e for the invariant mass in the A region (W = 1.232 GeV) and for the fixed high momentum 
transfers. The box correction in the SPT case has a weak dependence on e and reaches its largest value of —0.5% for 
Q 2 = 7.0GeV 2 . On the contrary, the correction calculated with the help of the SPMT and FM approaches shows a 
very strong dependence on e, which is enhanced in the region of the higher momentum transfers. For all approaches 
considered here, the biggest absolute value of the correction is observed for the case of backward electron scattering. 
If the invariant mass is above the mass of A resonance for Q 2 = 3.0 GeV 2 ~ 7.0 GeV 2 , we observe that the absolute 
value of the correction tends to become smaller in the case of backward electron scattering for the SPMT and FM 
approaches only. We find that the correction calculated using the SPT prescription is not sensitive to the changes of 
the invariant mass in the region of 1.232 GeV ~ 1.72 GeV. For the low momentum transfer (see Fig.© (right plot)), 
as in the case of tt° electroproduction, the absolute value of the correction obtained in the SPMT and FM approaches 
decreases substantially, hence justifying its strong Q 2 dependence. For the low momentum transfer of Q 2 = 0.4 GeV 2 , 
the dependencies on the invariant mass W, polar and azimuthal angles of the detected ir + are shown in Fig.©. In 
contrast with the high momentum transfers, for Q 2 = 0.4 GeV 2 the correction becomes rather small on the absolute 
scale and all the approaches agree with each other within 6%. We find that the correction in Fig.© (left plot) with 
fixed low momentum transfer and Ei a b = 1.645 GeV, is more sensitive on the relative scale to the variations of the 
invariant mass in the SPMT and FM approaches compared to the case of high momentum transfer. For the overall 
comparison between corrections obtained at the high and low fixed momentum transfers, we show dependencies of 
the two-photon box correction on the hadronic angular variables in Fig. © and Fig. QlUp , respectively. 

From Fig.© and (|10|1 it is evident that the box radiative correction has a visible effect on the angular distribution 
in the SPMT and FM approaches only. 

Another implementation of the box radiative corrections can be found in the studies of the tt + elastic electric form 
factor. Once again, in order to extract the Born cross section from the experimental data, it is essential to apply 
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Figure 7: Two-photon box correction (for detected 7r + ) dependencies on virtual photon degree of polarization parameter e for 
Q 2 = 3.0 GeV 2 (left plot), Q 2 = 7.0 GeV 2 (middle plot) and Q 2 = 0.4 GeV 2 (right plot). All plots are given for (j> 4 = 90° and 
9 4 = 90° and W = 1.232 GeV. Dot-dashed curve - SPT, dotted curve - SPT with qtt subtracted, dashed curve - SPMT, solid 
curve - FM approach. 
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Figure 8: ir + electroproduction box correction for fixed Q 2 — OAGeV 2 and Ei a t, = 1.645 GeV. Left plot: dependence on W 
with (j>4 — 90° and 84, = 90°. Middle plot: dependence on cos 64 with <f>4 = 90° and W = 1.232 GeV . Right plot: dependence 
on $4 with 6 — 90° and W = 1.232 GeV . Dot-dashed curve - SPT, dotted curve - SPT with air subtracted, dashed curve - 
SPMT, solid curve - FM approach. 




Figure 9: Examples of the two-photon box correction angular dependencies for all three approaches at fixed high momentum 
transfer, invariant mass and electron lab energy. All plots are given for Q 2 = 6.0 GeV 2 , W = 3.2 GeV and E iab = 9.2 GeV. 
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Figure 10: Examples of the two-photon box correction angular dependencies for all three approaches at fixed low momentum 
transfer, invariant mass and electron lab energy. All plots are given for Q 2 = 1.7 GeV 2 , W = 1.4 GeV and E iah = 5.754 GeV. 
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Figure 11: tt + electroproduction azimuthal (left plot, 64 — > 0°) and virtual photon degree of polarization (right plot, (f>4 = 90°and 
64 — s- 0°) dependencies of the two photon box correction for the fixed Q 2 = 6.0 GeV 2 , Ei a b = 10.9 GeV (for the case of azimuthal 
dependence) and W = 3.2 GeV. 

the radiative corrections to the measured cross section. For the highest momentum transfer of Q 2 — 6.0 GeV 2 and 
the invariant mass W = 3.2 GeV relevant to the proposed Ftt [5] experiment, Fig. lfTTj) shows the box correction for 
the case of forward tt + electroproduction (#4 — > 0°) in the center of mass of final hadrons. For all approaches, the 
azimuthal dependence (see Fig. (ITU (left plot)) of the correction is unsubstantial and has a variation of the order of 
2% over the range of 4>i € [0, 2vr]. The range of the correction is in the interval of —18% ~ —2.5% and the largest 
absolute value is observed for the correction calculated in the SPMT approach. Dependencies of the box correction on 
the degree of polarization parameter e are demonstrated in the Fig. fTTj) (right plot). It is quite interesting to see that 
the correction obtained in the SPMT approach exhibits very strong dependence to the variations of e for the case of 
7r + production in the forward direction. Overall, the correction calculated in the SPMT approach shows the largest 
sensitivity to e in the forward direction of the produced charged pion if compared to the SPT and FM approaches. The 
fact that the correction calculated using the FM method does not grow substantially, even for the case of backward 
scattering (e — ^0), for the case of deep-virtual kinematics (Q 2 = 6.0 GeV 2 , Ei a b — 10.9 GeV and W — 3.2 GeV) is 
most likely evidence of the strong suppression of the two-photon box amplitude by the monopole form factor in this 
case. The more detailed numerical results for the two-photon box correction relevant to the planned Fir experiment 
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are given in the Appendix. 



VI. CONCLUSION AND DISCUSSION 

In this work, we have evaluated the two-photon exchange correction to pion electroproduction. Calculations were 
performed with and without the soft-photon approximation. Two versions of the soft-photon approximation were used. 
In the first version, we applied the soft-photon approximation to the box amplitude consistently in the numerator 
and denominator algebra, hence preserving the hadronic current and satisfying the Ward identity. In the second 
version, we have followed the proposal of Maximon and Tjon (2l| and applied the soft-photon approximation to 
the numerator algebra only. To test the validity of the soft-photon approximation, we also calculated the two- 
photon box amplitude exactly, with no approximations, but instead we had to rely on a model for the reaction 
dynamics, namely to include the monopole form factor in all hadronic couplings. We found that the approximation 
Re [Co Tni\ 1 {—kj, JTij})] » —Co({ki,mi},{kj,nij}) used in [20] for the three-point scalar integrals resulted in a 
systematic shift in the value of correction by —2.3%. After comparing all three approaches, it became evident that the 
SPT method has produced a very weak kinematic dependence of the correction. Hence, the soft-photon approximation 
in the first approach removes the full kinematic dependence from the correction because the momentum of one of 
the virtual photons is neglected in the fermion denominators. This effect is mostly noticeable at higher momentum 
transfers. On the contrary, the SPMT and FM approaches do conserve the full kinematic dependence which is 
enhanced at higher momentum transfers. We found that the soft-photon exchange approximation provides a reasonable 
estimate for the self-consistent, model-independent calculation, but contributions from hard-photon exchange need 
to be included. It is also noted that various recipes for separating the soft-photon exchange contributions may lead 
to rather different numerical results. One may avoid the use of soft-photon approximation at a cost of introducing 
models for the two-photon exchange amplitude. Such a model (Formfactor-Model) was used in this work for illustrative 
purposes. It is realized that more sophisticated models are required in order to describe the pion electroproduction 
in the wide range of kinematics considered here. 

In conclusion, we found that the two-photon exchange contributions to the pion electroproduction may play a 
significant role for high momentum transfers and need to be included in the data analysis of precision experiments in 
electron scattering. 
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VII. APPENDIX 



Numerical results for the two-photon box correction for the planned Fn experiment are demonstrated in the table 
TblUfor the 6» 4 -> 0°. 
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Table I: Two-photon 7r + electroproduction box correction for the proposed Fir experiment with the 4>4 = 90° and 84 — > 0°. 



[11] I. Akushevich, N. Shumeiko, A. Soroko, Eur. Phys. J. C 10, 681 (1999). 
[12] I. Akushevich, Eur. Phys. J. C 8, 457 (1999). 

[13] A. Afanasev, I. Akushevich, V. Burkert and K. Joo, Phys. Rev. D 66, 074004 (2002). 

[14] D. Y. Bardin and N. M. Shumeiko, Nucl. Phys. B 127, 242 (1977). 

[15] P. G. Blunden, W. Melnitchouk, J. A. Tjon, Phys. Rev. C bf 81, 018202 (2010). 

[16] P.G. Blunden,W. Melnitchouk, and J. A. Tjon, Phys. Rev. C 72, 034612 (2005). 

[17] Y.-B. Dong, S.D. Wang, Phys. Lett. B 684, 123 (2010). 

[18] D. Borisyuk, A. Kobushkin, Phys. Rev. C 83, 025203 (2011). 

[19] T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118, 153 (1999). 

[20] Y.S. Tsai, Phys. Rev. 122, 1898 (1961). 

[21] L.C. Maximon, J.A. Tjon, Phys.Rev. C 62, 054320 (2000). 

[22] E. A. Kurae v et al, Phys. Rev. D 7 4, 013003 (2006). 

[23] A. Manohar, |a"rXiv:hep-ph/9305298| (1993). 

[24] A. Afansev, I. Akushevich, V. Burkert, K. Joo, Phys. Rev. D 66, (2002) 074004 

[25] A. N. Villano et al., Phys. Rev. C80 (2009) 035203 

[26] K. Joo et al. (CLAS), Phys. Rev. Lett. 88, 122001, (2002). 

[27] K. Joo et al. (CLAS), Phys. Rev. C 68, 032201 (2003). 

[28] A. Biselli et al. (CLAS), Phys. Rev. C 68, 035201 (2003). 



18 

[29] K. Joo ct al. (CLAS), Phys. Rev. C 70, 042201 (2004). 



